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We use QCD sum rules to estimate the leading-order universal form factors describ- 
. ing the semileptonic B decay into orbital excited D-wave charmed doublets, including 

O ■ the 2~) states {D^, D'^) and the (2~, 3~) states (1)2) D^)- The decay rates we 

predict are Ts^Dieu = ^B^D'.ei, = 2.4 x lO-i^GeV, TB^D2eu = 6.2 x lO-^^GeV, 



D ; and TB^D*eu = 8.6 x lO^^'^GeV. The branching ratios are B{B Dliu) = B{B 

; D'^ep) = 6.0 X 10-^ B{B D2^) = 1.5 X 10~^ and B{B D^iU) = 2.1 x 10~^ 

CN . respectively. 

^ ■ PACS numbers: 14.40.-n, 11.55.Hx, 12.38.Lg, 12.39.Hg 
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I. INTRODUCTION 

Higher excitations than D^*^ play an important role in the understanding of semileptonic 
B decays. Knowledge of these processes is important to reduce the uncertainties of the 



I measurements on other semileptonic B decays, and thus the determination of the Cabibbo- 

ic decay 
ation of 



^ I Kobayashi-Maskawa matrix elements, such as |V"cf,|. Theoretically, the semileptonic decay 

0\ '. processes are described by some form factors. The challenge for theory is the calcu 
O 



these decay form factors. Fortunately, the heavy quark effective theory (HQET) l|, with 
an expansion in terms of l/mq for hadrons containing a single heavy quark, provides a 
systematic method for investigating such processes. In HQET the approximate symmetries 
allow one to organize the spectrum of heavy mesons according to parity P and total angular 
momentum si of the light degree of freedom. Coupling the spin of the light degrees of 
freedom si with the spin of a heavy quark sq = 1/2 yields a doublet of meson states with 
a total spin s = s/ ± 1/2. For charmed mesons, the lowest lying states (0~, 1~) doublet 
{D, D*) are S-wave states with the spin of light degrees si = 1/2. The P-wave excitation 
corresponds to two series of states, one is the = 1/2 series, the (0"^, 1+) doublet (-Dq, D[); 
the other is the si = 3/2 series, the (1+, 2+) doublet {Di, D"^)- For D-wave states, those are 
(1~,2~) and (2^,3") doublets ((-DJ, D'2) and {D2, -D3)), corresponding to the spin of light 
degrees of freedom si = 3/2 and si = 5/2. The early study of the heavy-light mesons can be 
found in Ref. [2l. The S'-wave and P-wave charmed states have been observed so far. The 
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properties of these states have been extensively studied using different approaches during 
the past few years, including masses jl, 0], decay constants 0, 0, 0], and decay widths 
jsl, 0, l3, 111- For the Z)-wave charmed mesons, their properties were investigated with the 
potential model and QCD sum rules 12]. 

Semileptonic B decay into an excited heavy meson has been observed in experiments 
fisl . Q- Recently, BABAR has measured semileptonic B decays into orbitally excited 
charmed mesons Di(2420) and 1^2(2460) 15|. They also reported two new Dg states 



Dsj{2860) and Dsj{2690) in the DK channel, which may fit in the D-wave charm-strange 



doublets [16]. A similar state Dsj{2715) has also been observed by Belle [l7j. It is expected 
that the nonstrange D-wave charmed mesons will be found, and the measurements of the 
semileptonic B decays into these states become available in the near future. To this end 
we study the predictions of HQET for semileptonic B decays to D-wave charmed mesons. 

The semileptonic decay rate of a B meson transition into an charmed meson is deter- 
mined by the corresponding matrix elements of the weak axial- vector and vector currents. 
In the heavy quark limit these elements are described, respectively, by one universal Isgur- 



Wise function at the leading order of heavy quark expansion 18|| . The universal Isgur- Wise 



function is a nonperturbtive parameter. It must be calculated in some nonperturbative 



approaches. The main theoretical approaches are QCD sum rules [l9l], constituent quark 
models, and lattice QCD. The inve stig ations of semileptonic B decays into charmed mesons 
can be found in Refs. j^, [isl, 20, 21, 22] with different methods. In this work, we esti- 
mate the leading-order Isgur- Wise functions describing the decays B {D\,D'2)tP and 
B — > Dl)iV and give a prediction for the widths of the decays. 

The remainder of this paper is organized as follows. In Sec. [TT]we present the formulas of 
weak current matrix elements and decay rates. In Sec. IIIII we give the relevant sum rules 
for two-point correlators, and then deduce the three-point sum rules for the Isgur- Wise 
functions. Section [IV] is devoted to numerical results and discussions. 

II. ANALYTIC FORMULATIONS FOR SEMILEPTONIC DECAY 
AMPLITUDES B {Dl,D'^)iV AND B {D2,D^)il7 

The heavy-light meson doublets can be expressed conveniently by effective operators 



231 ]. For the ground doublet, the operator is 



Ha = ^[D;Y-D^,]. (1) 

The effective operators describing the meson doublets 1^(1^,2^) and D{2^,3^) are given 
by 

= ^-^[DTl^lu - DlJlig^'^ - + vn)], (2) 
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and 

H''' = ^-^[DT"!. - ^Jll.Df{gi!,g^, - ^g'',{Y - v^) - ^gi^iY - vn)]. (3) 

In these operators, D*, D, D^^, D*^^^'' , and 1^2^ separately represent annihilation 

operators of the Qq mesons with appropriate quantum numbers and |^ = f -7, f is the heavy 
meson velocity. The theoretical description of semileptonic decays involves the matrix 
elements of vector and axial-vector currents {V^ = c'-f'^b and = 07^756) between B 
mesons and excited D mesons. For the processes B —>■ [Dl, D2)£v and B — > {D2, Dl)iV, 
;hese matrix elements can be parametrized through applying the trace formalism as follows 



23|: 



{Dl{v\e)\{V-Ar\B{v)) = J^,yrfmiB^niy)[e* ■ viv'^ - ^^v'^ 



CvM, (4) 



{D'2{v,e)\{V - Ar\B{v)) = ^/mfn5^T,{yy>'W''' [v " 1) " v'v'^ + le^^'^^v'^vpl (5) 
and 



+ ^^-^^e^'^^v,v^], (6) 

{d;{v',6)\{v - Ar\Biv)) = ,/mm^MyKf^x^''^^[9'\^ + y)- + ^^"''^^p^II (7) 

where {V ~ A)'^ = C7'^(l — 75)6 is the weak current, y = v ■ v and Ti(y), T2{y) are the 
universal form factors, and e* , e*^, e'^px the polarization tensors of these mesons. The 
differential decay rates are calculated by making use of the formulas (jl]) to ([7]) given above: 
dT G%ylm'i.m%., . 3 

— (i?^Dt£77) = ^^-^(ri(y)r(y-l).(y + l).[(l + r?)(2y+l)-2ri(y2 + y+l)], 
^ ^ (8) 

"-{B D',tV) = 727r3 H nivmy -l)Hy + l)Hil + rlK4y -1) -2r2{3y' -y + 1)], 

^{B D^eu) = ^^^30!k^r2{ymy-l)H 

^ (10) 
dT G%y'im%,m'i., ^ . , 

— {B DltU) = 360^3 H T2{ymy-l)Hy+l)n{l+rl){lly+3)-2n{8y'+3y+3)], 

^ mo ^^^^ 

with Tj = {Di = Dl, D2, D2, for i = 1,2,3,4 ). In the equations above, we 

have presented the decay rates of B semileptonic decay processes B — > {Dl,D'2)tn and 

B — >• {D2, Dl)tP in terms of the universal form factors Ti{y) and T2{y), respectively. The 

only unknown factors in these equations are Ti{y) and T2{y), which need to be determined 

by nonperturbative methods. 
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III. SUM RULES FOR ISGUR-WISE FUNCTIONS 

In the calculation of Isgur-Wise functions in HQET by means of QCD sum rule, the 
interpolating currents are potentially important. In Ref. 0], two series of interpolating 
currents with nice propertties were proposed: 

jf:^:r^ = K{x)T^^p\D^Mi^ (12) 

or 

--"^ = 7^.(x)r];;;-^^(^.J(-^)P..?(^) (13) 

where i = 1,2 corresponding to two series of doublets of the spin-parity (j + 

X)(-i)^+^] + l)^""*^-*^], respectively. Dt^ = — v^{v ■ D) is the transverse 

component of the covariant derivative with respect to the velocity of the meson and 

p{ai...a,}p^j = symmetrize{r°i-"^(L)^J - -gt^'^^g^, ,H4«3-a,} ^^^^^^ 

3 '^l"2 

with the transverse metric g'^^ = g""^ — For the doublets of spin-parity , (j + 

X^C-i)^"*"^] and [j + the expressions for F"^'""-' (D^^J have been explicitly given 

in \m as 



|±i7^(-^)^D- ■ ■ • D^UD-f - 2^7^ P.J, for J^"^)^"^ 

j^.irHYD^f . .-m, for (j + i)(-i)^'"^ 



^^5(_,),^a.^., . . . D:r\ for (j + 1)(-^)^ 

^^i-^)W-^f . . . D^^, (D- - ^7r^P.J, for J^~'^' 

where = 7^ — i^v^ is the transverse component of 7^ with respect to the heavy quark 
velocity. 

For the D-wave meson doublets with = | and = | , where j = I and j = 2, the 
currents are given by the following expressions: 



1 

"^2, -,3/2 ~ 

and 



4,-^3/2 = -^^T'^^'^'Kl^^t^.Dt.q, (16) 



(17) 



V2 

which correspond to Eq. (|T2|) . and corresponding to Eq. f[T3|) are 



4-,3/2 = - \ltPt)Ptq, (19) 
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4^3/2 = -^T'^^^^'Kl.lt.Dt^Ptq. (20) 



and 



(21) 



1 



4S/2 = -^T-^''^'^"/^.7t/.A.A.(-0Ag, (22) 

where is the generic velocity-dependent heavy quark effective field in HQET and q 
denotes the light quark field. The tensors T"^'^'^ and T°'^^'^'^'^ are used to symmetrize 
indices and are given by ^] 

T-''^'' = lig^g?'' + 9Ti') - (23) 



/ a/3 ut' Act i a/3 M"" Ai/ i a/3 i/o- A/i , aA M'^ /3o' i aA 

- T7{9t 9t 9t +9t 9t 9t +9t 9t 9t + 9t 9t 9t + 9t 9t 9t 
15 

+ 9?' 9^9?" + 9^' 9^9?" + 9?'9r9r + 9^t'9r97)- (24) 



Usually the currents with derivatives of the lowest order (1121) are used in the QCD sum 
rule approach. However, currents with derivatives of one order higher (fT3l) are also used 
in some conditions because in the nonrelativistic quark model there is a corresponding 
relation between the orbital angular momenta and the orders of derivatives in the space 
wave functions. As for the orbital D-wave mesons, which corresponding to derivatives of 
order two, it is reasonable to use the currents ( fTTl) . ( fTSl) . ( fT9l) and ( l20l) . 

These currents have nice properties, they have nonvanishing projection only to the 
corresponding states of the HQET in the mq oo limit, without mixing with states 
of the same quantum number but different si. Thus we can define one-particle-current 
couplings as follows: 

= r : {Dl{v,e)\r\Q) = h^/mBle*^, (25) 

= 2- : {D',{v,e)\rP\Q) = f^^^e*"^ (26) 

= 2- : {D^{v,e)\r^\Q) = /av^^*"^ (27) 

= 3- : (D*(t;,e)|J"^» = h^^e*^^\ (28) 

The couplings are low-energy parameters which are determined by the dynamics of 
the light degree of freedom. Since the pairs (/i, /g) and (/2, /s) are related by the spin 
symmetry, we will consider fi and /2 hereafter. The decay constants /j can be estimated 
from two-point sum rules, therefore we list the sum rules after the Borel transformation. 
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For the ground-state heavy mesons, the sum rule for the correlator of two heavy-light 
currents is well known. It is 20 1 



^-4^ -le^X .V"/^..-i(,-,)(i-9L), (29) 

For the sf = | doublet, when the currents (fT9l) and (!20|) are used, the corresponding sum 
rule is : 

/! s_e-'~^"'^^^ = r u'e-l^du - r u'e-Z^dui^GG). (30) 

'2 2%2 3x28^0 

For the sf = | doublet, when the currents (fT71) and (fT8|) are used, the corresponding sum 
rule is : 

•^-'5 5 X 2^TT^ Jo 3 X 26 Jo ^ vr ^ ^ ^ 

As we have just mentioned, for the amplitudes of the semileptonic decays into excited 
states in the infinite mass limit, the only unknown quantities in ([8]), ([9]), (fTOll and (fTT!) are 
the universal functions Ti{y) and T2{y). In Ref. |2J] the form factors ti(|/) and T2{y) were 
estimated through QCD sum rule by using currents with derivatives of lower order, ([T^ to 
( fT8|) . Considering that the corresponding relation between the orbital angular momentum 
and the order of the derivative mentioned above, we use the currents f|T9|) and fl20l) instead 
of f|T5|) and f|T6|) for the (D^, doublet. As for the {D2, -D3) doublet, we also use the 
currents (IT7|) and (ITSI) . 

In order to calculate this two form factors by QCD sum rules, we study the analytic 
properties of three-point correlators: 



rfWze'('='-^-'=-^)(0|T[Ji^(x)J{;5'''')(0)JoV(^)|0) = T{uj,J,y)C^yy, (32) 
d'xd'ze^^''--'-\mjt^^)4T'\^)4A'^)\^) = r'(a;,u;',i/)£tr^, (33) 



where Jy^'^ = h{v')'y^h{v) and J'^^^'^ ^ = h{v')'^^'^^h{v). The variables k{= P — rrihv) 
and k' {= P' — rricv') denote residual "off-shell" momenta of the initial and final meson 
states, respectively. For heavy quarks in bound states they are typically of order Aqcd 
and remain finite in the heavy quark limit. T{uj, cu' , y) and T'ioj, uj\ y) are analytic functions 
in the "off-shell" energies u = 2v ■ k and u' = 2v' ■ k' with discontinuities for positive values 
of these variables. They also depend on the velocity transfer y = v ■ v', which is fixed in a 
physical region. Cv,a are Lorentz structures. 

Following the standard QCD sum rule procedure, the calculations of T{uj,u;',y) and 
T'{u,uj\y) are straightforward. First, we saturate Eqs. flS^ and fl55]) with physical inter- 
mediate states in HQET and find that the hadronic representations of the correlators as 
follows: 

/ ' N f-,if-jMy) 

^hadron{uJ, ,y) = --j . ; + higher resonances, (34) 

(2A_ 1 - a; - te){2A_,j^ - u - le) 
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where f-j^ are the decay constants defined in Eqs. (l25il and ( l27j) . A__j; = — uiq. 

Second, the functions can be approximated by a perturbative calculation supplemented 
by nonperturbative power corrections proportional to the vacuum condensates which are 
treated as phenomenological parameters. The perturbative contribution can be repre- 
sented by a double dispersion integral in u and u' plus possible subtraction terms. So the 
theoretical expression for the correlator has the form 

dudiy , / , ' — ; + subtractions + r'^°"'^(a;, , w). (35) 

[u — uj — tejlu — uj — le) 

The perturbative part of the spectral density can be calculated straightforward. Confining 
us to the leading order of perturbation, the perturbative spectral densities of the two sum 
rules for Tiiy) and T2{y) are 

p,,rM y\v) = 24?-— ^^^'[(3^^ - (1 + 2y)i2uu' - u'')] 

(y-i)2 

xQ{u)Q{u)Q{2yuu' -u"^ -u"^), (36) 

and 

Ppertii^, v\ y) = ^-—T-^ ^[(5z/ - Uyu' + 3z/>2 + (3z/ + z/')(22/2 - 2y + 1^^] 

xe{u)e{u')e{2yuu' -u"" -u'^). (37) 
Following the arguments in Refs. [s], 25], the perturbative and the hadronic spectral densi- 



ties cannot be locally dual to each other, the necessary way to restore duality is to integrate 
the spectral densities over the "off-diagonal" variable z/_ = z/ — z/' , keeping the "diagonal" 
variable = ^^i^ fixed. It is in z/+ that the quark-hadron duality is assumed for the 
integrated spectral densities. The integration region can be expressed in terms of the vari- 
ables h'_ and and we choose the triangular region defined by the bounds: < z/+ < Uc, 
—2sJ~^^h'+ < y- < 2^1^ As discussed in Refs. js, [isl, the upper limit uj^ for z/+ in 

the region |[(?/ + 1) — \/ y"^ — l]i^co ^ ^ |(^^co + ^c2) is reasonable. A double Borel trans- 
formation in UJ and uj' is performed on both sides of the sum rules, in which for simplicity 



we take the Borel parameters equal j^. 



24j : Ti = T2 = 2T. In the calculation, we have 



considered the operators of dimension D < 5 in OPE. After adding the nonperturbative 
parts, we obtain the sum rules for Ti and T2 as follows: 

n(y)/-,V2/-,3/.e-(--/.^A.3/.)/T = dv,e-^vX 

^ + (38) 



3 X 25 (y + 1)2 ^ TT 



r2(l/)/-,i/2/-,5/2e-(^->^/^+'^->^/^)/^ = ^2lVT7Xi \ ^^^^ 



!i 4 



7r2 (1 + yY J, + + 
^ ^ {^GG). (39) 



3 X 24 (y + 1)3 ^ TT 



8 



We also derive the sum rule for T2 by using the currents (12T|) and (122|) . which appears to 
be 



%-25 g ° 
+ 3x24 (y + l)3<V^^^- 

IV. NUMERICAL RESULTS AND DISCUSSIONS 

We now evaluate the sum rules numerically. For the QCD parameters entering the 
theoretical expressions, we take the standard values: (gg) = — (0.24)^GeV^, (a^GG) = 



0.04GeV'^, and ml = O.SGeV^. In the numerical calculations, we take 2.83GeV [21, Il0| for 
the mass of the si = 5/2 doublet and 2.78GeV for the si = 3/2 doublet. For mass of initial 
B meson, we use tub = 5.279GeV [2Q ] . 

In order to obtain information of Ti{y) and T2{y) with less systematic uncertainties in 
the calculation, we divide the three-point sum rules by the square roots of relevant two- 
point sum rules, as many authors did [s], [20, 24 1, to reduce the number of input parameters 
and improve stabilities. Then we obtain expressions for the Ti{y) and T2{y) as functions 
of the Borel parameter T and the continuum thresholds. Imposing usual criteria for the 
upper and lower bounds of the Borel parameter, we found they have a common sum rule 
"window": 0.7GeV < T < l.SGeV, which overlaps with those of two-point sum rules (l29il . 
f l30|) and fl3T|) (see Fig. 1). Notice that the Borel parameter in the sum rules for three-point 
correlators is twice the Borel parameter in the sum rules for the two-point correlators. In 
the evaluation we have taken 2.0GeV < ujco < 2.4GeV 0,0, 2.8GeV < ujci < 3.2GeV, and 
3.2GeV < uJc2 < 3.6GeV. The regions of these continuum thresholds are fixed by analyzing 
the corresponding two-point sum rules. According to the discussion in Sec. Illlt we can fix 
oj'^ and ojc in the regions 2.3GeV < uj'^< 2.6GeV and 2.5GeV < Uc < 2.7GeV. The results 
are showed in Fig. 2. The resulting curves for Ti{y) and T2{y) can be parametrized by 
the linear approximation 

Ti{y)=T,{l)[l-pl^{y-l)l ri(l) = 0.14 ±0.03, p^^ = 0.13 ± 0.02; (41) 

r2{y)=T2{l)[l-pl^{y-l)l r2(l) = 0.57 ±0.09, p^^ = 0.78 ± 0.13. (42) 

The errors mainly come from the uncertainty due to uo^s and T. It is difficult to estimate 
these systematic errors which are brought in by the quark-hadron duality. The maximal 
values of y are yrAx = ymlx = (1 ± 2)/2'^i,2 ~ 1-213 and y^l^ = ymlx = (1 ± r| 4)/2r3,4 ^ 
1.201. By using the parameters Vet = 0.04, Gp = 1.166 x 10~^GeV~^, we get the semilep- 
';onic decay rates of B {D\,D'2)tn and B {D2, D^)il7. Consider that tb = 1.638ps 



26| , we get the branching ratios, respectively. All these results are listed in Table [H 

Because of the large background from B — »• D'^*^tu decays, there is no experimental data 
available so far. As we can see from Table HI the rates of semileptonic B decay into the 
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FIG. 1: Dependence of Ti{y) and T2{y) on Borel parameter T at y = 1. 




FIG. 2: Prediction for the Isgur-Wise functions Ti{y) and T2{y). 



sf = I doublet are tiny and our results are larger than those predicted by Ref. [24] in the 
B to s[ = I charmed doublet channels. The difference comes because the way in which 
we choose the parameters is different from theirs. They chose the parameters according 
to other theoretical approaches. In contrast, we choose the parameters following the way 
of Ref. j^. In addition, we also estimate the universal form factor T2{y) with the sum 
rule (HOj) and we get almost the same result as fj42|) . When trying to estimate the Ti{y) by 
using the currents (HM and (ITBI) . we find that after the quark-hadron duality are assumed 
the integral over the perturbative spectral density becomes zero. As for the P-wave and 
the F-wave mesons, similar results can be obtained after the calculations above have been 
carefully repeated. 

The semileptonic and leptonic B decay rate is about 10.9% of the total B decay rate, 
in which the S'-wave charmed mesons D and D* contribute about 8.65% 26|] and the 
P-wave charmed mesons contribute about 0.9% [20]. Our results then suggest that the 
D-wave charmed mesons contribute about 0.04% of the total B decay rate. Sum up the 
branching ratios of these semileptonic B decay processes, the eight lightest charmed mesons 
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TABLE I: Predictions for the decay widths and branching ratios 



Decay mode Decay width T (GeV) Branching ratio Branching ratio of Ref . [241 
B D\tV 2.4 X 10"^^ 6.0 x 10"^ 

B D'^iU 2.4 X 10"^^ 6.0 x 10"^ 

B D2tV 6.2 X lO"^'^ 1.5 X lO"'' 1 x 10"^ 

B D*£u 8.6 X lO"^'^ 2.1 x lO""^ 1 x 10"^ 



contribute about 9.59% of the B decay rate. Therefore, semileptonic decays into higher 
excited states and nonresonant multibody channels should be about 1.31% of the B decay 
rate. Whatsoever, our result is just a leading-order estimate of the contribution of the 
D-wave charmed mesons channels to the semileptonic B decay. 

In summary, we estimate the leading-order universal form factors describing the B meson 
of ground-state transition into orbital excited D-wave charmed resonances, the (1~, 2~) 
states {Dl, D'2), which belong to the sf = | heavy quark doublet and the (2^, 3^) states 
{D2, -D3), which belong to the sf = | heavy quark doublet, by use of QCD sum rules 
within the framework of HQET. The semileptonic decay widths as well as the branching 
ratios we get are shown in Table [11 The predictions are larger than those predicted by Ref. 



2J] . This needs future experiments for clarification. We also prove that when sf = | the 



interpolating currents (|T2|) and f|T3|) proposed in Ref. [j] are really equivalent. It is worth 
noting that in the estimate of the semileptonic B decay form factors when the currents 
( IT2l) with quantum numbers of light degree of freedom s[ = | , |^ are used for the 
excited charmed mesons, we find the perturbative contributions vanish after the quark- 
hadron duality are assumed. In this case we should use the currents (fT3l) which contain 
derivatives of one order higher. 
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